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Abstract 

We show that anisotropic negative effective dispersion relation can be achieved in pure dielectric 
rod-type metamaterials by turning from the symmetry of a square lattice to that of a rectangu- 
lar one, i.e. by breaking the rotation symmetry of effective homogeneous medium. Theoretical 
predictions and conclusions are verified by both numerical calculations and computer based sim- 
ulations. The proposed anisotropic metamaterial, is used to construct a refocusing slab-lens and 
a subdiffraction hyperlens. The all-dielectric origin makes it more straightforward to address loss 
and scaling, two major issues of metallic structures, thus facilitating future applications in both 
the terahertz and optical range. 
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In recent years, metamaterials (MMs) made of a stack of small structural elements have 
attracted substantial attention. The inherent ability of high-index structures to support 
strong subwavelength resonances allows the associated MMs to behave as effective media 
with various unusual macroscopic electromagnetic (EM) properties, such as negative per- 
mittivity/permeability and negative refraction lsj]. With peculiar collective EM charac- 
teristics, MMs can be used to fabricate innovative devices like subwavelength refocusing 
lenses |4[ and invisibility cloaks pi]. 

Currently, one of the most important issues on MMs is to scale the MM unitcells down in 
size to allow operating a yet higher frequencies, i.e. terahertz and optical regions js, 7|. In 
this context, dielectrics rather than metals might be preferred by MMs, since the dispersive 
metals would be non-ideal and exhibit significant loss, thus compromising MMs design and 
application at high frequencies. Till now, great efforts have been devoted to developing 
homogenization theory for the design and fabrication of pure dielectric MMs. Rod-type 
dielectric resonators with high permittivity have been used to fabricate left-handed (LH) 
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MMs in the microwave range and the terahertz range [8|, |9| . Recent scaling analysis have 
also pointed out that LH behaviors at optical frequencies would also be possible by use 
of rod-type MMs made of silicon [10]. Rod-type dielectric MMs with an isotropic density 
of rods (supported e.g. by an ordered or disordered square-lattice of rods) will commonly 
possess an isotropic dispersion in the two dimensions jsl, Q. While isotropic dispersion 
is of advantageous in a number of cases, there is however also a strong call for realizing 
anisotropic MMs for practical applications, as is in general required for a transformation 
medium. 

In this paper, we report on both theoretical and numerical studies of the dielectric rod- 
type MMs with asymmetrical shaped unit cells. General closed-form expressions for the 
effective constitutive parameters are obtained from a macroscopic viewpoint, i.e. in the 
long wavelength limit. We show that the dispersion relationship of such dielectric MMs 
would be naturally spatially anisotropic with negative refraction, supported by parameters 
extracted from numerical computation and EM phenomena in other simulations. By using 
this kind of anisotropic rod-type MMs, a farfield slab-lens with different refocusing distance 
is conveniently realized. Moreover, this kind of MMs can also be implemented to fabricate 
a slab-type hyperlens with sub-diffraction limit. Both of these designs are supported by 
full-wave simulations. 
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We start by considering the array composed of standard cylindrical resonators, as shown 
in Fig. Ufa). For simplicity, we emphasize the ordered case, but the composite could be 
disordered, since the low-frequency macroscopic EM properties of this kind metamaterial do 
not rely on band structure phenomena, as they are mainly determined by the low-lying Mie 
resonances of the isolated cylinders with large permittivity 8, Q. As discussed in Ref. js], 
the displacement current distribution inside the cylinder in the basic two Mie resonance 
modes would be similar to that of the conductive current in metallic rod and Split-ring 
Resonator (SRR), thus also supporting effective negative permittivity/permeability of the 
metamaterial. Fig. QJb) shows the electric field and magnetic field distributions associated 
with two first basic Mie resonance modes, respectively. When the dielectric cylinders are 
arranged, in a regular lattice or in a disordered fashion, the whole composite will exhibit 
isotropic EM responses provided that the average density of rods is the same for e.g. the x 
and the y directions. Spatial anisotropy may be introduced by relaxing this isotropic density, 
thus breaking the effective rotational symmetry of the long-wavelength effective medium. 

Here, without loss of generality, rod-type dielectric metamaterial would have a rectangular 
shaped unit cell, as shown in Fig. Q^c). Basically, to treat this kind of rod-type composite 
as an effective medium, the typical dimensions of a single unit cell have to be much smaller 
than the wavelength in free space [8j. The effective parameters, i.e. e e g and ^ eff , are 
defined in Ref. |8J as (X?) = e e g(E} and (i?) = ^ efr (H"). From Mie scattering theory, 
all the correlated fields can formally be expressed by superposition of vector cylindrical 
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111 ] . However, for the rods array in two dimensions, when the exciting electric field 



is polarized along the axis, i.e. the z-direction in Fig. HJc), the practical geometry of the 
lattice complicates a closed-form analytical solution for the integration over the rectangular 
shaped area S\. Although numerical evaluations can make precise estimations, we will for 
the ease of analysis approximate the integral over the rectangular domain by performing 
the integrals on an equivalent region formed by So, S[ and S' 2 , see Fig. QJd). In this way 
simplify the problem, yet still capturing the correct physics. With these approximations, if 
a is much smaller than the free wavelength in vacuum, in which case high-order scattering 
is very weak, straightforward algebraic operation would eventually give 
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As an example, we include numerical evaluations of the above expressions for a practi- 
cal case. We assume that the dielectric rods are made of Bao.sSro.sTiOs (BST) ceramic, 
which has a dielectric constant around 600 in the microwave range [§J. It has previously 
been shown that a MM made from rods of this ceramic will have a pass band with negative 
refractive index 0, For the dimensions of the asymmetric lattice we consider Z=7.4mm 
and h=7mm, which will support an anisotropic dispersion relation, as compared to isotropic 
property of the corresponding square lattice. Fig. |2fa) shows the calculated effective re- 
fractive index n c s = ■^/e^g-Jp^g in both x and y directions. Considering the pre-conditions 
of our homogenization, we only focus on those frequencies at which the absolute value of 
the refractive index is small, see the inset of Fig. |2]^a). The difference between n x and n y 
indicates the anisotropic property of this MMs. In addition, we have also performed wave 
simulations and extract the effective refractive index in the two directions by a retrieval 
method [12]. The extracted data is shown in Fig. [2(b). Obviously, Fig. [2ta) is qualitatively 
consistent with Fig. [21(b). 

For an effective medium, the dispersion relation of an effective medium could formally be 
written in an analytical form like ^f- + ^f- = k$, with n x and n y being the refractive index 
in x and y directions, respectively. Though the different effective refractive index associated 
with different directions is responsible for the anisotropic EM response, it is quite difficult to 
verify the complete dispersion curve by a retrieval method as proposed by Ref. 12] because 
of the rotational asymmetry. Here, we instead illustrate the anisotropic property, by the 
example of a device which require an anisotropic MM, i.e. refocusing of a point source by a 
slab lens. 

As known, a slab lens with refractive index of -1 can make sub-wavelength refocusing in 



free space 4|. As an effect of a compensating bilayer, a general slab lens with sub- wavelength 
resolution can be realized by anisotropic negative refractive index metamaterials (NIMs) 
whose dispersion relation has a form of k% + % = k%, where r = (d' — d)/d is relating to the 
distance between source and its image (d') vs. width of the slab lens (d) [l3]. By using the 
rod-type dielectric metamaterial, such a slab lens would be easily designed and constructed. 
For simplicity, ordered metamaterials are adopted here. The parameters for the first slab 
lens are a=l mm, 1—7 mm, and h—7 mm, with the rod's relative dielectric constant of 600. 
Since the dispersion relation of this slab lens would be isotropic, a refocused image can be 
found at a distance d' = 2d, see the power flow density distribution (PFDD) in Fig. [3(a). 
Next, we keep I unchanged, but vary h in the y direction. In case (J), h is increased to 7.4 
mm, and from the PFDD in Fig. [3(c) we see that the source is refocused again at d' < 2d 
with a little down shifting of operation frequency. In case (II), h is decreased to 6.4 mm, 
in which case the refocused image is found at a distance d' > 2d. For comparison, we have 
also simulated the PFDD of refocusing by an ideal negative index anisotropic slab lens, see 
Fig. [3(b) (isotropic material with n ~ —1), and panels [3(d) and [3(f) is for an anisotropic 
material with k\ + \A 2 ky = k\ and k 2 x + = k%, respectively. 

Clearly, from Fig. [31 the refocusing of a MM slablens is governed by the intrinsic effective 
negative constitutive relations. However, it is very hard to obtain subwavelength resolution 
in such an imaging case, since its effective impedance is generally mismatched with that of 
the free space [8|, leading to inevitable reflection at the interface. Fortunately, there is an 
alternative way to achieve sub- wavelength resolution with the current MMs. From Fig. [5J 
the effective refractive index in one direction will turn imaginary at some frequencies just 
beyond the upper bound of the double negative band, i.e. n x is imaginary at 4.78 GHz 
in Fig. [2(b), while n y stays real. Thus, the MMs' effective dispersion would be hyperbolic 
in a narrow frequency band, e.g. from 4.77 GHz to 4.79 GHz in Fig. [2(b). With these 
characteristics, this MM can be used to form a hyperlens which can distinguish sources with 
subwavelength separation placed close to the interface 
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Again, BST cylinders with 1 mm radius are adopted to construct a metamaterial with 
1=7 mm and h=6A mm. In this way, a 0.128 m wide slab hyperlens is designed, see 
Fig. E(a). In the simulation, two point sources with a separation of 16 mm are placed in 
front of the hyperlens at y=-0.005 m. The simulated power flow density distribution around 
4.81 GHz is shown in Fig. E(a), and indeed we observe two images at the interface. For a 
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clearer visualization, panel H^b) shows cross-section cuts of the PFDD at y=-0.005 m and 
y=0.128 m. Keeping the radiation wavelength in mind, the ability to distinguish the two 
sources and the image indicates sub-diffraction imaging by the MMs hyperlens. 

In conclusion, we have proposed a way to design anisotropic negative-index MMs by 
means of pure dielectric rod-type rectangular lattices. Results are supported by a homoge- 
nization approach as well as by numerical calculation and simulations. From our theoretical 
analysis, the effective dispersion is represented by closed form expressions like a realistic 
anisotropic media existing in nature. 

With the anisotropic MMs proposed here, we have constructed a far field refocusing 
slablens and a sub diffract ion hyperlens. Since the effective constitutive parameters play the 
key role in achieving the EM phenomena, anisotropic rod-type MMs could also be used to 
some other applications, e.g. constructing an invisibility cloak or transformation medium. 
Moreover, our study in this paper points out flexible design of dielectric rod-type MMs, 
which avoids the loss and scaling issue of metallic structures in realizing negative properties, 
indicating many potential applications in terahertz and optical range in future. 
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(c) (d) 

FIG. 1: (color), (a) The configuration of rod-type dielectric metamaterials. (b) The fundamental 
electric and magnetic resonant modes in a dielectric cylinder with z polarization incidence; (c) The 
cross section schematic of rod-type metamaterial unit cell and (d) its simplified integral region. 
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FIG. 2: Effective refractive index in different directions of anisotropic rod-type metamaterial. 
(a) Theoretical estimation; (b)Data extracted from numerical simulation. In the calculation and 
simulation, a=lmm, !=7.4mm, /i=7mm, and the relative dielectric constant of the rod is 600. 
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FIG. 3: PFDD in the refocusing by three different slab lenses made of pure dielectric metamaterials. 
The periodicities in x and y directions are (a)7mmx7mm,(c)7mmx7.4mm,(e)7mmx6.4mm; the 
widths of the three superlens are (a)0. 126m, (c)0. 1258m, (e)0. 128m. In the three simulations, the 
operation frequencies change slightly, respectively. The insets of (a), (c) and (e) are the cross 
section of the power flow density distribution at y=0.21m. Subfigures (b), (d) and (f) are the 
PFDD in the refocusing by three ideal anisotropic slablens with negative refractive index. The 
region (0 < y < 0.128m) is occupied by the slablens for (b), (d) and (f). 
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FIG. 4: Sub-diffraction by a slab-shape hyperlens made of rod-type dielectric metamaterial. 
(a)PFDD; (b) Cross section draws of PFDD at y=-0.005m and y=0.128m. In the simulation, 
/ w4.81GHz, a=lmm, /=7mm, /i=6.4mm, while the interval of the two sources is 16mm. In (b), 
the data at y=0.128m is normalized by a factor of 100 to have a better comparison. 
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